A theoretical scheme of thermal-hght ghost imaging by A^th-order intensity correlation 
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In this paper, we propose a tlieoretical scheme of ghost imaging in terms of A''th-order correlated 
thermal light. We obtain the Gaussian thin lens equations in the ghost imaging protocol. We show 
that it is possible to produce N — 1 ghost images of an object at different places in a nonlocal fashion 
by means of a higher-order correlated imaging process with an TVth-order correlated thermal source 
and correlation measurement. We investigate the visibility of the ghost images in the scheme, and 
obtain the upper bounds of the visibility for the A^'th-order correlated thermal-light ghost imaging. 
It is found that the visibility of the ghost images can be dramatically enhanced when the order of 
correlation becomes larger. 
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I. INTRODUCTION 

Ghost imaging with thermal Hght Jl[lJlJl[i,[E0,[i, 
i, M, [m, Elli M, M, 113, m El la ISl has been 
studied extensively in recent years. Bennink and cowork- 
ers [l| first pointed out that ghost imaging can also be 
realized using a classical source with the appropriate cor- 
relations. A thermal or quasi-thermal source can exhibit 
such a classical correlation. A very close formal analogy 
was demonstrated between ghost imaging with thermal 
and quantum-entangled beams in Refs. [3, S S , which 
implied that classically correlated beams were able to 
emulate the relevant features of quantum ghost imaging. 
Number of experiments on thermal light ghost imaging 
have been performed in Refs. [1, [ifl, [ul, [13, [111. A 
unified treatment of classical and quantum ghost imag- 
ing was established in terms of Gaussian-state analysis 
inRef. [H. 

Recently, some attention has been paid to the physics 
of thermal-light ghost imaging [H, [13, [H, [H, [H [H, 
\2T\ . [H, [2^ and higher-order-coherence or correlation ef- 
fects of thermal light [10, [Ml, [13, [H [H, [H . The higher- 
order coherence or correlation has shown attractive prop- 
erties in practical applications. Multiphoton imaging 
with thermal light is one of these exciting areas. In a 
previous work, our group theoretically initiated a study 
of thermal-light ghost imaging in terms of higher-order 
correlated thermal light 30). Wc proposed a thermal- 
light ghost imaging scheme with third-order correlated 
thermal light. In this scheme, a third-order correlated 
thermal-light source, a test optical arm and two refer- 
ence optical arms are used to produce two ghost images. 
Two ghost images are created at two different places in a 
nonlocal fashion as a consequence of the third order cor- 
relation of the involved optical fields. It was shown that 
the third-order correlated imaging includes richer cor- 
related imaging effects than the second-order correlated 
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one. In this paper, we want to propose a thermal-light 
ghost imaging scheme with A'^th-order intensity correla- 
tion and investigate the visibility of the ghost images. 
We show that the usual second- and third-order ghost 
imaging are only two particular examples of our present 
scheme. 

The paper is organized as follows. In Sec. II, begin 
with the thermal-light source with A^th-order correlation, 
we will propose the thermal-light ghost imaging scheme 
with A^th-order intensity correlation. In Sec. Ill, the 
visibility of the ghost images in the proposed scheme will 
be investigated. Finally, wc shall conclude our paper with 
discussions and remarks in the last section. 



II. GHOST IMAGING SCHEME WITH 
iV-ORDER CORRELATED THERMAL LIGHT 

The basic setup for ghost imaging with A^-order cor- 
related thermal light is indicated in Figure 1 which in- 
cludes one test optical arm and A^ — 1 reference optical 
arms. First arm is the test arm. An unknown object with 
transmission function T[x) and a collective lens with fo- 
cal length fc are placed on the test arm while there is one 
imaging lens with focal lengths fi on each reference arm. 
The object is placed at the focal plane of the collective 
lens at distance zi from the thermal-light source. The 
distance between ith imaging lens and the thermal-light 
source is z^o with i = 2,3, ■ ■ ■ , N. The first detector Di 
is a bucket detector placed at the focal plane of the col- 
lective lens on the right-hand side. The other detectors 
are scanning detectors Di with i = 2,3, - ■ ■ ,N which are 
placed at distances zn from ith imaging lens. The sig- 
nals from the A^ photon counting detectors are sent to 
an electronic coincidence circuit to measure the rate of 
coincidence counts. 

The thermal light source with A^th-order intensity 
correlation, usually obtained by illuminating a laser 
beam into a slowly rotating ground glass, is divided 
into A^ beams, which can be implemented by an ap- 
propriate combination of A^ — 1 beam splitters. Con- 
sider a monochromatic plane wave described by the field 
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Eq exp[i(fcoz — Wot)] illuminating a material containing 
disordered scattering centers. After scattering, the field 
can be written as E{x, z,t) = J E{q) ex.p[i{q ■ x + kzZ — 
u'ot)](iq where q is the transverse wave vector intro- 
duced by the random scattering and satisfies the relation 
|qp + — fc§. Hence, E{q) is a stochastic variable obey- 
ing Gaussian statistics. However, the scattered waves 
with different transverse wave vectors are statistically in- 
dependent. If |q| ^ fco, the scattered field can be approx- 
imately written as E(x,z,t) = exp[(i(fcoZ — w^t)] 
where A{x) = / ^(q) exp[j(q- x)]dq is the slowly varying 
envelope. As a result, we have defined a monochromatic 
thermal light random in both strength and propagation 
direction. According to the Wiener-Khintchine theorem, 
the first-order spectral correlation satisfies the following 
expression 

(i?*(q)ii;(q')>=^(q)^(q-q'), (1) 

where S'(q) is the power spectrum of the spatial fre- 
quency. For any field with thermal statistics, all high- 
order correlations can be expressed in terms of the first- 
order ones due to {E{ci}) = [sSl. Then, A^th-order 
spectral correlation of thermal light can be written as 
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■^(qjv), (2) 



where the prime in the summation means that all of re- 
peating terms are subtracted from the summation over 
Ti and Vj . After passing through a combination of A^ — 1 
beam splitters, the thermal light with A^th-order spec- 
tral correlation is divided into A^ correlated thermal- 
light beams which are input light beams of the test and 
reference arms. For simplicity, we consider the one- 
dimensional case, and let xq and x„ with n = 1, 2, • • • , A^ 
be the transverse coordinates of the source plane and de- 
tection planes, respectively. Here we have assumed that 
the starting planes of the A^ thermal-light beams are over- 
lapped. Let li\{x\^xo) be the impulse response function 
of the test arm, and i?„(a;„,a:o) with n = 2,3,- ••,A^ 
be the impulse response function of nth reference arm, 
respectively. Assume that the light field on the nth de- 
tection planes is denoted as E{xn), which is connected 
with the optical field of the thermal-light source through 
the following relation 



(3) 



E{Xn)^ / hn{Xn,-qn)E{qn)dqn 



where /i„(a;„,(7„) = (l/\/27r) / a;o) exp(-zg„a;o)dxo 

is the Fourier transformation of the impulse response 
function _ff„(x„,a;o). Then, the A^th-order correlation 
function of the joint intensity at the A^ detection planes 
can be expressed as 
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Xi,X2, ■ ■ 



■,XN) = {l[E*{x,)Eix'r)). (4) 



The rate of coincidence counts is governed by the A^th- 
order correlation function given by equation ([4]) which 
can be calculated in terms of the impulse response func- 
tions of the relevant optical systems. The impulse re- 
sponse functions ^] in the reference arms can be written 
as 
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where Zri ^ fr with r — 2, 3, • • • , A^, and the impulse 
response function in the test arm is given by 
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FIG. 1: (Color online) Schematic of the setup for implement- 
ing ghost imaging with A^th-order correlated thermal light. It 
consists of one test arm and A — 1 reference arms. The test 
arm includes the object T{x), the collection lens Fc, and the 
bucket detector Di. The object and the detector are placed in 
its two focal planes of the collection lens, respectively. Each 
reference arm includes one imaging lens and a scanning detec- 
tor denoted by Fi and Di with i = 2, 3, ■ ■ ■ , A — 1, respectively. 

In order to obtain the A^th-order correlation function 
given by equation (4), we first calculate the following 
integrations: 



S {q)\hr ixr,-q)\ dq, 



(8) 



S {q)h*{xr,-q)hr' {xr',-q)dq, (9) 
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where both r and r' may take 1, 2, • • • , iV, but r ^ r' . In 
the broadband Hmit, S(q) can be regarded as a constant 
S(0). Hence we have the total intensity of the thermal 
light Sq = J S{q)dq ~ S{0)qa, where go is the spectral 
bandwidth of the source. Substituting the impulse re- 
sponse functions of the test and reference arms given by 
equations ^ and ([7]) into equations ([5]) and (O, we ob- 
tain 

r ^(0)fc f 2 , J frS{0)qo 

(11) 
(12) 



Crl 




27r V ifr - Zrl) ifr' ~ Zr'l) 



where both r and r' may take 2, 3, • • • , but r ^ r', and 
we have introduced two phase functions 



(Prl{x, q) = —\ZrO-Zi + 



kxi 



+qXr - k{ZrO + Zrl - - 2/c), 



(13) 



2 ifr — Zrl) 2 (fr' — ^r'l) 



/ ^ _^ ^ ^r'l/r Zrl fr 

-q {Xr' - Xr) + k{Zr>0 + Zr'l - ZrO - Zrl), 

(14) 

where the scaled positions Xr — frXr/{fr — Zri). 

From equations (11) and (13), it is straightforward 
to see that when the positions of the object, ghost im- 
ages, and the lenses obey the following Gaussian thin lens 
equations for the A^th-order correlated imaging 



1 1 



ZrO ^1 ^rl fr 



(r = 2,3,---,iV), (15) 



the cross-correlation functions of intensity fluctuation 
can be simplified as 



Crl — 



kfr 



m / 

2tT Y i27r/c {fr - Zrl) 



frfr' 



m / 
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expi-!^)T(X.), 



(16) 

(0) J(X^, 

(17) 



When Gaussian thin lens equations are satisfied, and 
^2 7^ ^3 ■ ■ ■ 7^ the A^th-order correlation function 
of the joint intensity at the N detection planes given by 



equation Q can be reduced to the following form 

G^^^ {xi,X2, ■ ■ ■ ,Xn) 
N 

= ^ (-^n-fra • • ■ ^rjv_2 |C'rjv„ilP) +hh---lN, 
ri^rj=2 

(18) 

where the prime on the right-hand side of above equation 
means that repeating terms in the summation over and 
rj should be subtracted. The last term on the right-hand 
equation (18) is the background term which is the multi- 
plication of the intensity distribution at the N detectors. 
It does not contribute to the correlated imaging, but it 
may affect the visibility of produced ghost images. Each 
term in the summation over and Tj is the multiplica- 
tion of the intensity distribution at detector Di and the 
intensity fluctuation correlation between Di and Di with 
i = 2, 3, • • • , A^, which gives the information of the object 
imaged at D^. 

In particular, when the A^ — 1 reference arms are iden- 
tic, i.e., /2 = /a = • • • = /at, Z20 = 230 = • • • = znq and 
2:21 = Z31 = • • ■ = ZMi, we have h = h = ■ ■ ■ ^ In = I , 
A^th-order correlation function becomes 



N 



(19) 



where we have introduced the parameter x = 

S^{0)kfr/[8n^fc{fr-Zrl)]. 

Equation (19) indicates that for the object placed at 
the test arm, A^ — 1 ghost images can be produced in 
the A^ — 1 reference arms through coincidence count 
measurements upon A^ detectors. Each ghost image 
is an amplified image of the object with the amplify- 
ing rate /,./ {fr — Zri) with r = 2, 3, • • • , A^ since Xr = 
frXr/{fr — Zrl). From the Gaussian thin lens equations of 
A^th-order correlated imaging given by equation (15) we 
can see that the higher-order correlated imaging exhibits 
richer imaging effects. In fact, the imaging equation (15) 
indicates that Zro — zi is the object distance while Zri is 
the image distance for rth joint path with r = 2, 3, • • • , A^. 
Just as the ordinary imaging law, when the object dis- 
tance of the rth joint path Zro — zi is greater (less) than 
the focal length fr, the correlated image is real (virtual). 
As a specific example, we consider the case of A^ = 4. In 
this case, from the imaging equation (15) we can obtain 
the following eight ghost-image configurations: (1) when 
ZrO — zi > fr with r = 2, 3, and 4, the three ghost images 
are real; (2) when Zro — zi < fr with r = 2, 3, and 4, the 
three ghost images are virtual; (3) when Z20 — zi > f^,, 
^30 — zi < fs and Z40 — zi < /4, one ghost image is real 
at the position X2, the other two ghost images are vir- 
tual at the positions X3 and X4; (4) when Z30 — zi > f^, 
Z20 — zi < fi and Z40 — zi < f4, the ghost image at the 
position X3, two ghost images at the positions X2 and X4 
are virtual; (5) when Z40 — zi > fi, Z20 — zi < /2 and 
^30 — zi < fa, the ghost image at the position X4 is real. 
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two ghost images at the positions X2 and ^3 are virtual; 
(6) when Z20 - > /2, ^30 - < /a and Z40 - zi < f^, 
two ghost images at the positions X2 and 2:3 are real, 
the ghost image at the position is virtual; (7) when 
Z20 - zi> f2, Z4Q - zi < /4 and Z30 - zi < fs, two ghost 
images at the positions X2 and x^ are real, the ghost im- 
age at the position 2:3 is virtual; (8) when Z30 — zi > fa, 
Z40 — zi < and Z20 — zi < /2, two ghost images at 
the positions 0:3 and X4 are real, the ghost image at the 
position X2 is virtual. 

It is obvious that when = 2, 3, the usual results of 
the second- and third-order thermal-light ghost imaging 
are recovered. 



III. VISIBILITY OF GHOST IMAGES 

The visibility is an important parameter to estimate 
the quality of a ghost image. We first generalize the 
visibility of the second-order correlated imaging 
to the case of A^th-order correlated imaging as follows 



V 



(N) 



(N) 



(20) 



For our present scheme of the A^th order thermal-light 
imaging, substituting equations (19) into equation (20), 
we obtain the visibility 



El2\TiX^)\' 



qoj\Tix)\Ux+j:l2\nx^)\ 



(21) 



which indicates that the visibility depends on not only 
the spectral bandwidth of the thermal-light source qo and 
the transmission function T{x) but also the order num- 
ber of the correlation of the thermal-light source. The 
visibility can be enhanced with the decrease of the spec- 
tral bandwidth. An increase of the transmission leads to 
a decrease of the visibility since when the transmission 
area increases, more points contribute to the background 
that directly makes the visibility decrease. 

In order to see the relationship between the visibility 
and the order number of the correlation, making use of 
equations (10) and (11), we have 



{h){Ir) 



\{E*{x,)E{Xr)W 



^ qo \T {x)\^ dx, (22) 



where r — 2,3, ■ ■ ■ , N . According to the Cauchy- 
Schwartz inequality (/i) (1^) > \{E* (xi) E {xr))\^, which 
means that go / \T {x)\^ dx > 1. Hence, we find that 







max 
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N-1 



N 



(23) 



which implies that the visibility of the ghost images can 
be dramatically enhanced when the order of correlation 



becomes larger. For A'th order correlated thermal-light 
imaging, the upper bounds of the visibility is given by 
V^^^ = {N~ 1)/N. As expected, when A^ = 2 the upper 
bound of the visibility of the second correlated imaging 
M is 1/2- 



IV. CONCLUDING REMARKS 

In conclusion, we have proposed a theoretical scheme 
of ghost imaging in terms of A^th-order correlated ther- 
mal light. Our scheme includes one test arm and A" — 1 
references arms. One object is placed on the test arm. 
N —1 ghost images are produced in the reference arms in 
a nonlocal fashion by means of a higher-order correlated 
imaging process with an A^th-order correlated thermal 
source and correlation measurements. We have derived 
the Gaussian thin lens equations which the positions of 
the ghost images obey in the ghost imaging protocol. 
We have also investigated the visibility of the ghost im- 
ages in the scheme, and obtained the upper bounds of 
the visibility for the A^th-order correlated thermal-light 
ghost imaging. It has been shown that the visibility de- 
pends on not only the spectral bandwidth of the thermal- 
light source and the transmission area of the object but 
also the order number of the correlation of the thermal- 
light source. It is found that the visibility of the ghost 
imagines can be dramatically enhanced when the order 
of correlation becomes larger. Our present scheme is a 
many-ghost imaging protocol. On one hand, it gives rise 
to a theoretical origin for developing many-ghost imag- 
ing technology. This gives rise to the possibility of ex- 
perimentally producing correlated many ghost images. 
In fact, the higher-order correlated imaging opens up 
new avenues for realizing multi-port information process- 
ing. On the other hand, physically these ghost images 
stem from higher-order coherence or correlation of opti- 
cal fields. In this sense, the appearance of the many ghost 
images reveals an observable physical effect of higher- 
order coherence or correlation of optical fields. Hence, it 
is of very significance to study higher-order correlated 
imaging not only for well understanding the essential 
physics behind the higher coherence or correlation of op- 
tical fields but also for developing multi-port information 
processing technology. The experimental realization for 
the many ghost imaging protocol proposed here and prac- 
tical applications of the many ghost imaging phenomenon 
deserves further investigation. 
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